Technical Newsletter No. 8 



APPLIED SCIENCE DIVISION 



IBM 



APPLIED SCIENCE DIVISION 

Technical Newsletter No, 8 

September 1954 



This issue of the Newsletter is devoted solely to articles on the 
IBM Type 650 Magnetic Drum Data Processing Machine. Three 
of the articles were written by George R. Trimble, Jr. , and 
Elmer C. Kubie who were members of the IBM Mathematical 
Planning Group at Endicott, N.Y. Mr. Trimble is also the 
author of the remaining three articles. All the routines and 
programs presented in this issue have been tested and checked 
out on the Type 650 in Endicott. 



Copyright, 1954, by International Business Machines Corporation 
590 Madison Avenue, New York 22, New York 



CONTENTS 



1. Principles of Optimum Programming the IBM Type 650. . . .5 
G. R. Trimble, Jr. and E. C. Kubie 



2. An Interpretative Floating Decimal System for the IBM 

Type 650 1"^ 

G. R. Trimble, Jr. and E. C. Kubie 



3. Floating Decimal Sub-Routines for the IBM Type 650 ... . 37 
G. R. Trimble, Jr. 



4. IBM Type 650 Loading Routines 44 

G. R. Trimble, Jr. and E. C. Kubie 



5. A Method for Performing Double Precision Arithmetic 

on the IBM Type 650 -60 

G. R. Trimble, Jr. 



6. A Method for Performing Complex Arithmetic on the 

IBM Type 650 67 

G. R. Trimble, Jr. 



PRINCIPLES OF OPTIMUM PROGRAMMING 
THE IBM TYPE 650 

G. R. Trimble, Jr. 
E. C. Kubie 



Introduction 

The IBM Type 650 Is a parallel- serial calculator utilizing a magnetic drum for 
storage. It has been designed with "ease of use" as one of the primary considerations. 
The programmer is not burdened with timing restrictions which he must always keep in 
mind. Interlocks make it impossible to violate timing conditions in such a way as to 
cause the machine to give erroneous results. If, however, proper recognition and 
analysis of the sequence of events occurring within the machine is made, one can signifi- 
cantly increase the overall speed by properly locating data and instructions. Optimum 
programming is the technique by which data and instructions are located in such a manner 
as to minimize or eliminate, if possible, non-productive waiting or seeching time. 

The basic cycle of the 650 is the "word time" or time required to read one word. 
As there are 50 words written around the drum, each word time is equal to 1/50 of a 
drum revolution. Since the drum revolves at 12, 500 revolutions per minute a word 
time is equal to . 096 milliseconds. Each operation can be analyzed in terms of word 
times. The object here is to analyze each operation to see how many word times are 
required in its interpretation and execution. In this manner the basic or fundamental 
word times are determined and from these a set of rules is derived through which 
optimum programming may be effected. 

Angular drum locations are all that must be considered when optimum coding. 
The 650 is completely synchronous so that the words located in equivalent angular 
positions are in phase and are equivalent from a timing viewpoint. Thus, the words 
in addresses n modulo 50 are equivalent, so that, for example, addresses 0003, 
0053, 0103, . . . , 1903, 1953 are in phase and therefore, equivalent so far as 
optimum programming is concerned. This means that there are 40 different 
locations along the length of the drum, any of which may be used without loss 
of time. 

The lower half of the accumulator can be read into or out of only during an 
even word time and the upper half of the accumulator only during an odd word time. 
Thus, when the operation called for is one which uses the accumulator, it may be 
necessary to wait for an even word time or an odd word time in order to obtain the 



data. This is why the cycles 'wait for even', or 'wait for odd' on the sequence chart, 
are required. 

The D-I address system used in the 650 greatly facilitates optimum programming. 
Data and instructions may be in any locations on the drum and instructions may be 
taken in any desired sequence. As mentioned above, once the optimum location has 
been determined there are 40 different locations in which the datum or instruction 
could be placed since the 40 bands are in phase. 

It is also significant that the address of many instructions will be optimum 
within a range of angular locations. This is due to the 650's ability to overlap 
arithmetic execution with search for the next instruction. 

The characteristics of the 650 which facilitate optimum programming are, 

1. It is synchronous, that is, all timing is controlled from timing pulses 
on the drum. 

2. Interlocks make it impossible to cause errors by violation of timing 
conditions. 

3. Equivalence of angular locations in the 40 bands. 

4. Parallel operation, that is, simultaneous execution of operation and 
search for next instruction. 

5. The D-I address system which makes possible flexible location of data 
and instructions. 

Sequence Chart 

The sequence chart shows the steps taken in the interpretation and execution 
of each operation. It is drawn in segments which are usually equivalent to one 
word time. Where events are shown in parallel, they are performed simultaneously. 
Most operations branch into parallel execution shortly after the data is made available 
to the arithmetic unit. One branch indicates the arithmetic process and the other 
the obtaining of the next instruction. 

Analyzing the steps carried out by the machine in the performance of any 
operation, one sees that in every case certain fundamental word times are required. 
Begin consideration of each instruction at the time when the instruction has been 
located but has not entered the program register. Starting at this point, the first 
word time of every operation is used to transfer the instruction from its memory 
location to the program register. The next word time is used to initiate the 
interpretation of this instruction. This is done by transferring the data address to 
the address register and the operation code to the operation register. The steps 
performed during the word times beyond this point will depend upon the particular 
instruction being executed. 

Consider the shift right 3 operation. The portion of the sequence chart for 
this particular operation is shown in Figure 1. Assume that the shift instruction 
is in location n and that n is odd. During word time n, the instruction is read 



into the program register. During word time n+1 the operation code (30) is 
transferred to the operation register and the data address (0003) is transferred 
to the address register. Word time n+2 "Enable Shift Control", is used to set 
up the necessary control circuits to perform the shifting. Since n+3 is even, 
the next cycle, "Wait for Even, " is skipped and shifting begins immediately. 
At this point parallel operation begins and the actual shifting process on the 
right branch occurs simultaneously with the process of searching for and obtaining 
the next instruction indicated on the left branch. 

The restart signal which occurs during word time n+3 indicates that the 
data address in the address register is no longer needed. During word time 
n+4 the instruction address is transferred to the address register. Word time 
n+5 is the "Enable Program Register Read In" cycle which indicates that a new 
instruction is to be read into the program register. Thus, beginning with word 
time n+6 the 650 starts searching for the next instruction. This search may 
require from word times to as many as 49 word times depending on the location 
of the next instruction. 
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FIGURE 1 



Examining the arithmetic branch it is seen that during word times n+3 throuijh 
n+8 the actual shifting takes place. The interlock at point A is removed during 
word time n+9. The function of this interlock will be explained later. 

Although the arithmetic portion of the operation is not completed until word 
time n+9, the program register is able to accept the next instruction at word 
time n+6. K the next instruction were placed in a location corresponding to 
word times n+6, n+7, or n+8, it would be read into the program register before 
completion of the arithmetic portion of the shift operation. In addition to this, 
interpretation of this instruction would begin in that the operation code would be 
transferred to the operation register and the data address transferred to the 
address register. Thus, the first two cycles of the interpretation of the next 
instruction occur in parallel with the execution of the shift instruction. At 
this point, interpretation must cease since further steps may make use of 
control circuits or portions of the arithmetic unit which are already in use 
(unless it is a read or punch instruction). The interlock at A is provided to 
insure that interpretation of the next instruction does not proceed past 
this point. 

If the next instruction were placed in a location corresponding to word 
time n+9, the instruction would be read into the program register during 
completion of the arithmetic portion of the shift operation. In this case only 
the first cycle of the interpretation of the next instruction occurs in parallel 
with the execution of the shift instruction. 

If the next instruction were placed in a location corresponding to word 
time n+10 it would be read into the program register immediately upon 
completion of the shift operation. 

Each of the locations corresponding to word times n+6 through n+10 
effectively reduces the search time to since there is no wait time between 
completion of the shifting operation and the obtaining of the next instruction. 
However, if the next instruction were placed in a location corresponding to 
word times n+6 through n+9, one cycle of the interpretation of the next instruction 
would take place, if it were restricted to a location corresponding to word times 
n+6 through n+8, two cycles of interpretation of the next instruction would take 
place. I 

Thus, not only can the search for the next instruction take place during 
completion of the arithmetic portion of an operation, .but also interpretation 
of this next instruction can begin. If the next instruction calls for a read or punch 
operation, execution of either of these operations can proceed since they do not 
require use of the arithmetic unit. In this case the interlock A does not stop 
interpretation of the next instruction after the first 2 cycles but it will be noted 
that read and punch instructions bypass this interlock. 



As a numerical example assume n=643, that is, the instruction is in location 
0643. The Instruction address should be 0649, 0650, or 0651 to reduce the access 
time to zero. Of course, locations 0699, 0700, 0701, etc. , are equally good. 

There are three classes of operations. They are read, punch and arithmetic 
operations. Correspondingly, there are three interlocks. The read interlock "R" 
stops execution of the following read instructions, until completion of the previous 
read instruction. If the following instructions are punch or arithmetic operations 
however, they may proceed without delay. The punch interlock "P" stops execution 
of the following punch instructions, until completion of the previous punch instruction. 
If the following instructions are read or arithmetic operations however, they may 
proceed without delay. The arithmetic interlock "A" stops execution of the following 
arithmetic instructions until completion of the previous arithmetic operation. If the 
following instructions are read or punch instructions however, they may proceed 
without delay. Thus simultaneous reading, punching, and computing is possible. 
However, simultaneous execution of two arithmetic operations, two punching 
operations or two reading operations is not possible. 

Referring to the shift right three example again, it is seen that since read or 
punch instructions are not held up by the interlock at A, they should be placed in a 
location corresponding to word time n+6 (649 in the numerical example). If the 
next instruction is another arithmetic operation, there is no advantage to n+6 or 
n+7 over n+8, since the interlock at A causes the machine to wait for completion of 
the shifting process. 

The rules developed for optimum coding are based on an effective step being 
executed during every word time. This assumption is made to keep the rules 
simple. For this reason, it is best that the next instruction be placed in a location 
corresponding to word time n+8 (651 in the numerical example) if it is another 
arithmetic operation. This allows maximum overlap and assures an effective step 
for every word time (no waiting at interlock points) . 

Because of the above considerations the rules for shift instructions are stated 
in terms of 2 limits. The lower limit is best if the next instruction calls for a read 
or punch operation. The upper limit is best if the next instruction is an arithmetic 
operation. If the upper limit is not convenient, any location within.the two limits will 
result in zero access time. However, when considering the succeeding instructions 
and their data or instruction addresses, one should apply the rules as if the upper 
limit were used. With this in mind, the following rules apply for the shift right 3 
instruction. 



Shift Right Three 


Location of Next Instruction | 


Liower Limit 


Upper Limit 


n even 


i=n+7 


i=n+9 


n odd 


i=n+6 


i=n+8 



This example demonstrates how the sequence chart is used to determine optimum 



locations. It would be too lengthy a task to show how every rule is obtained from the 
sequence chart, so we will state the rules and leave the proofs for those readers 
sufficiently curious. 



Rules for Optimum Programming 

n : location of instruction 

d : data address of instruction 

i : instruction address of instruction 

U '• lower limit of i 



Sm : sum of multiplier digits 
Sq : sum of quotient digits 



^u- 



upper limit of i 

location of argument found by table look up operation 



Operation 

Add, subtract, etc., 10,11,15,16, 

17,18,60,61,65,66,67,68 
Multiply, 19 
Divide, 14^64 

Store Lower Accumulator, 20 
Store Upper Accumulator, 21 
Load Dist. , Store Dist. , 69, 24 
Store D Address, Store I Address, 22, 


n even 
d=n+ 


n odd 
d=n+ 


d even 
i=d+ 


d odd 
i=d+ 


3 
3 
3 
5 
4 
3 
23 3 


3 
3 
3 
4 
5 
3 
4 


5 
21+2Z;m 
61+2Sq 

3 

3 

3 

3 


4 
20+2Sm 
60+2Sq 

3 

3 

3 

3 




Operation 

Table Look Up, 84 


n even 
d=n+ 


n odd 
d=n+ 


a even 
i=a+ 


a odd 

i=a+ 


3 


3 


5 


6 



Operation 



44 



Branch Non-Zero Upper, 

Branch Non-Zero, 45 

Branch Minus, 46 

Branch Overflow, Branch Dist. 8, 47 

91-98 
Branch Dist. 8, 90, 99 
No Operation, Stop, 00, 01 



n even 
d=n+ 



"3" 
4 
3 

3 

4 



n odd 
d=n+ 



3 
3 

3 

4 



n even 
i=n+ 



4^ 

5 

4 

5 
5 

4 



n odd 
i=n+ 



~5" 
4 

4 

5 
5 

4 



No. of Positions 



"oar 



Shift Operations 

Shift Right, 30 
Shift Left, 35 
Shift and Count, 36 



r 



< 



0) 

(2) 
(3) 
(4) 
5) 
(6) 



V(9) 



Shift and Count Only, 36 (10) 



Shift and Round, 31 



(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
.(10) 



n even 



11 =n+ 



Jufn+ 



^^ 

7 

7 

9 
U 
13 
15 
17 
19 
21 

23 

7 

9 

11 
13 
15 
17 
19 
21 
23 
25 



il=n+ 
~5 

6 

6 

6 

6 

6 

6 

6 

6 

6 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 



ii,=n+ 



6 
6 
8 

10 
12 
14 
16 
18 
20 

22 

6 
8 

10 
12 
14 
16 
18 
20 
22 
24 
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800X Addresses 

Since the storage entry switches (8000), distributor (8001), lower accumulator 
(8002) and upper accumulator (8003) are addressable, special consideration must be 
given to the cases where one of these locations is addressed. The storage entry 
switches and distributor are always immediately accessible. However, the lower 
accumulator can be read into or out of only during an even word time. The upper 
accumulator can be read into or out of only during an odd word time. Thus, for 
purposes of optimum codiug the storage entry switches and the distributor may be 
treated as being equivalent to any address, the lower accumulator as equivalent 
to an even address and the upper accumulator as equivalent to an odd address. 

K, for example, a reset add upper operation (or any add tjnpe operation) has a 
data address of 8000 or 8001, these addresses can be treated as being equivalent 
to n+3 and the instruction address of that instruction can be determined accordingly. 
If a data address of 8002 were used, and the instruction were in an even location, 
the optimum location for the datum for that instruction would be Ji+3, which is odd 
in this case. Therefore, an extra cycle must be taken to wait for an even location 
so that the lower accumulator may be read out. Since the effective data address 
is then even, the rule i=d+5 must be used to determine the location of the next 
instruction. Similar analyses may be made for each of the other cases. 

The following table gives the rules for determining the instruction address of an 
instruction in cases in which a data address of 800X is used. 



Operation 


Data 
Address 


n even 
i=n+ 


n odd 
i=n+ 


Add, Subtract, etc., 
(10,11,15,16,17,18, 
60,61,65,66,67,68) 


8000 
8001 
8002 
8003 


7 
7 
9 
7 


00 CO 00 CO 


Load Distributor, 69 


8000 
8001 
8002 
8003 


6 
6 
7 
6 


6 
6 
6 

7 



Since the accumulator may be in use when an instruction address of 8002 or 8003 
is given, it would seem to be possible to take the next instruction from the accumulator 
before the arithmetic operation was completed. For example, if an 8002 instruction 
address were used on a multiply operation, one of the partial products might be taken 
as the next instruction rather than the final product. For this reason, an additional 
interlock is provided to prevent the next instruction from being taken from an 800X 
address until completion of the arithmetic portion of the operation. This prevents 
the above from occurring. 

The following table gives the rules for determlng the effective instruction address 
in those cases in which an instruction address of 800X is used. 
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Operation 


Inst. 
Address 


d even 
i^+ 


d odd 


Add, Subtract, etc. 

(10,11,15,16,17,18, 

60,61,65,66,67,68) 

No Complement Cycle Required 


8000 
8001 
8002 
8003 


5 
5 
6 

5 


4 
4 
5 
4 


Add, Subtract, etc., 
(10,11,15,16,17,18, 
60,61,65,66,67,68) 
Complement Cycle Required 


8000 
8001 
8002 
8003 


7 
7 
8 
7 


6 
6 
7 
6 


Load Distributor, 69 


8000 
8001 
8002 
8003 


3 
3 

4 
3 


3 
3 
3 

4 


•V means "equivalent to" 









By application of the principles used in determining the above rules it is possible 
to determine equivalent addresses for any other desired case. It is only necessary to 
remember that an 8002 address is equivalent to an even drum address and 8003 is 
equivalent to an odd drum address. 

Consider the instruction 1006438000. Since 0643 is odd, the rule in d+4 must be 
used. Therefore, the instruction at 8000 can be programmed as though it were located 
in location 0647. 
Techniques for Using Optimum Programming 

The gains obtained by optimum programming will depend upon the particular 
problem considered and the skill of the programmer. Every instruction cannot be 
optimum as conflict will occasionally exist between instructions. For example, 
data placed optimumly for a store operation may not be optimumly located for a later 
add operation. Furthermore, an instruction which is preceded by several branch 
instructions cannot usually be optimumly located for each of these branch instructions. 

Techniques for efficient use of optimum programming wiU be developed throi^h 
experience. It is too lengthy a task to thoroughly explore techniques for its use in 
this paper. However, certain general principles can be applied to any program one 
may wish to consider. Some of these principles will be indicated in the following 
paragraphs. 

The first consideration that must be made is when should optimum programming 
be used. 

Any problem in which the speed of input or output is appreciably reduced due to 
lengthy calculations, can justifiably be programmed optimumly. 

Once it has been determined that optimum programming is necessary, there are 
two ways in which it can be done. One could simply program the problem optimumly 
by straightforward application of the rules derived above without consideration of the 



tz 



program as a whole. Such a procedure may result in data and instructions being 
poorly located for the latter part of the problem. This type of programming can 
be done quickly and with very little more difficulty than if the instructions were 
programmed In sequence. Even though it is done roughly, it will usually result 
in a significant increase In the overall speed, and will be well worth the small 
amount of additional effort required. 

The second type of optimum programming is where the programming is done 
elegantly. This type of programming requires a greater amount of work and thought 
on the part of the programmer, instead of simply programming each step optimumly 
as it occurs, the programmer must think ahead to see how this might possibly 
affect later steps which wiU use the same data or perhaps branch to the same 
instruction. Many possibilities exist when programming in this manner. One can 
see after he has completed the problem how a simple re -arrangement of the 
beginning may improve the latter part of the problem. It may be necessary to re- 
program the same problem several times in order to obtain the most efficient 
program. Obviously such a procedure would require much more time than "sequential" 
programming or the "rough" optimum programming described above. 

Elegant optimum programming should be used only when it is desirable to use 
the 650 in the most efficient manner possible. The types of problems for which this 
is necessary are those which must be done over and over on a mass basis and which 
realize reduced input-output speeds. It will not always be necessary to program the 
entire problem this way but only those segments or sub-routines which are most 
frequently used. For example, sub-routines, such as floating decimal operations, 
square root, sine, and cosine e'valuations, which occur frequently in technical 
applications, should be programmed to function in the most efficient manner possible. 
Similar applications which will occur in commercial problems are extension from 
gross to net in payroll calculation, insurance dividend calculation, and bill computation 
in public utility customer accounting. 
Example 

The following example of programming for square root is included to help clarify 
how optimum programming is used. When programmed sequentially, this routine 
requires approximately .321 seconds (assuming 5 iti&rations). Secondly is shown 
the same routine coded optimumly. The time required is reduced to .152 seconds. 
Thus, the increase in speed is a factor of 2.1 to 1. Note, however, that the 
optimumly programmed routine requires 2 additional storage locations so that the 
constant 1/2 wUl always be located optimumly. Experience to date has shown that 
the gain realized may be a factor of less than 2 to 1 to a factor as large as 6 or 7 to 1. 
Computing During Input and Output Operations 

At a card punching rate of 100 cards per minute up to 544 milliseconds are available 
for computing. This is approximately 5600 word times or 110 drum revolutions. At 
a card reading speed of 200 cards per minute up to 25,7''-milliseconds are available for 
computing. (This is approximately 2700 word times or 54 drum revolutions.) 
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AN INTERPRETATIVE FLOATING DECIMAL SYSTEM FOR THE IBM TYPE 650 

G. R. Trimble, Jr. 
E. C. Kubie 

Introduction 

This floating decimal system will perform 18 basic operations using a floating 
decimal number system by means of interpretative programming. It was designed 
with coding convenience as a prime objective. This is illustrated by the fact that one 
of the basic operations is essentially a vector by vector multiplication. 

Floating decimal instructions, that is, instructions which are to be interpreted, 
have a negative sign. A floating decimal instruction consists of a two digit operation 
code, a 4 digit address specifying the location of the first factor, a 4 digit address 
specifying the location of the second factor (if required) and another 4 digit address 
specifying the location of the third factor (if required). It is a variable address system 
in that only as many addresses as are needed for the particular operation are required. 
Thus, some operations use only one address, some require two addresses and others 
require 3 addresses. 

Floating decimal instructions will be taken from consecutive memory locations. 
If an instruction requires one or two addresses, the instruction is stored in one memory 
location. If three addresses are required, the third address is stored in the memory 
location immediately following the one containing the instruction. 

Should a positive instruction appear, it is interpreted as a normal 650 instruction 
and subsequent instructions are not interpreted. Thus, upon occurrence of a positive 
instruction the 650 operates in its usual D-I mode. This continues until an instruction 
address of 0026 is given which causes control to return to the interpretative routine. 
The program will return to the floating decimal mode of operation at the point of 
departure. 

For example, consider the following sequence: 
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Location 
n 

n+1 

n+2 

n+3 



Contents 
floating point instruction (-) 

floating point instruction (-) 

floating point instruction (-) 

normal 650 instruction (+) 



(The contents of n+3 are interpreted as a normal 650 instruction including the instruction 
address for sequencing. Normal execution of instructions continues until 0026 is used 
as an instruction address, at which time the program returns to the floating decimal 
mode of operation beginning with the instruction in location n+4). 



Number Form 

The number form is as follows: 



XX' xxxxxxxx-t-^ 



-Exponent +50 
-Mantissa 



Instruction Form 

The instruction form is as follows: 



-floating point command (01-18) 
-1st address 
-2nd address 



location i 'xx "xxxx*'ixxxji-^ 
location i+1 oo xxxx oooo- 

i+1 has this form only if the third address is needed, otherwise, the form of i+1 is 

the same as the form of i. 



Floating Decimal Accumulator 

A floating point accumulator referred to as K is used for accumulation and 
accumulative multiplication. It does not need to be addressed in operations which make 
use of it. For example, operation 01 (A+B— -*K) uses two addresses, A and B. The 
sum is automatically stored in K as the result of this operation. 
Operations 

Below is a list of the operations Airtth their codes, addresses required and estimated 
average time of execution. 

Estimated Average Time 



Code and Addresses 

01, A, B 

02, A 

03, A, B 

04, A 

05, A, B, C 



Operation 
A+B ^K 



A+K- 



-^K 



A-B- 
K-A- 
A+B- 



->K 






62. 


8 


54. 


7 


63. 


1 


59 


5 


81 


8 
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Code and Addresses 

06, A, B, C 

07, A, B, C 

08, A, B 

09, A, B, C 

10, A, B 

11, A, B 

12, A 

13, A, B 

14, A, B 

15, A, B 

16, A, B 

17, A, B, C 

18, Ai, Bi, n 



Operation 
A-B-^ — VC 

AxB hC 

AxB >K 

A/B yC 

A/B ►K 

BRMIN A 
BR 

BRNZ A 
(AxB)+K-*K 
K-(AxB)->K 

vir — J.B 

V A'^+B^'+C >K 

n 

S A.B.-^K 
i=l ^ ^ 



Estimated Average Time 

96.2 

84.8 

78.4 

92.2 

72.2 

28.4 

18.7 

30.8 
102. 5 
102.5 
157.9 
416.3 

(32. 4+92. 3n) 



Explanation of programs 

1. General Interpretation 

The general interpretation routine takes instructions from consecutive memory 
locations and analyzes them to see if they are normal 650 instructions or if they must be 
interpreted. If the instruction is a normal 650 instruction, it is executed as such. If 
the instruction is to be interpreted, the routine obtains the factor at address A and stores 
it in location 0037. Control is then transferred to the proper sub-routine. The constant 
in location 0193 facilitates use of the translating routine TRl. The amount of translation 
is placed in the instruction address positions of this constant. This amount is then added 
to the operation code and control is transferred to the translated sub-routine. 

2. Addition Sub-Routine 

The addition sub-routine adds the factors in locations 0037, in which A is 
stored, and 0057, which is the floating decimal accumulator K. The result is normally 
stored in the floating decimal accumulator K. 

3. Multiplication Sub-Routine 

The multiplication sub-routine multiplies the factors in 0037, which is A, and 
0089, in which B is placed. The product is normally stored in 0037. 

4. Division Sub-Routine 

The division sub-routine divides the factor in 0037 which is A, by the factor 
B, which had previously been placed in the lower accumulator. The result is normally 
stored in 0057 which is the floating decimal accumulator. 
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The following sub-routines interpret the second and third addresses if required, 
obtain the necessary factors and modify the addition, multiplication and division sub- 
routines so as to store the result in the desired location. In some cases, they also 
modify portions of other sub- interpretative routines. 

5. Interpretation of 01 

This sub-interpretative routine obtains the factor B and transfers control to 
the addition sub-routine. 

6. Interpretation of 02 

Since the factors are already in place all that is necessary for this sub- 
interpretative routine is to transfer control to the addition sub-routine. 

7. Interpretation of 03 

This sub-interpretative routine obtains the factor B, reverses its sign and 
transfers control to the addition sub-routine. 

8. Interpretation of 04 

The sign of the factor A is reversed and control is transferred to the 
addition sub-routine. 

9. Interpretation of 05 

Since the addition sub-routine makes use of the floating decimal accumulator 
K, the contents of K must be temporarily transferred to another location, and returned 
after the operation has been completed. This sub-interpretative routine obtains B, 
modifies the addition sub-routine to store the result in location C and transfers control 
to the addition sub-routine. After completion of the operation, control is returned to the 
sub-taterpretative routine and the addition sub-routine is restored to normal. 

10. Interpretation of 06 

This sub-interpretative routine modifies the sub-interpretative routine, 
"Interpretation of 05" so that the sign of factor B is reversed. 

11. Interpretation of 07 

This sub-interpretative routine obtains the factor B, modifies the multiply 
sub-routine to store the product in location C and transfers control to the multiply 
sub-routine. After completion of the multiplication, control is returned to this sub- 
interpretative routine and the multiply sub-routine is restored to normal. 
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12 Interpretation of 08 

The factor B is obtained, the last instruction of the multiply sub-routine is 
modified to store the product in K and control is transferred to the multiply sub-routine. 
After completion of the multiplication, control is returned to the sub-interpretative 
routine "Interpretation of 07" which restores the last instruction of the multiply sub- 
routine to normal. 

13. Interpretation of 09 

The factor B is obtained, the last instruction of the divide sub-routine is 
modified to store the quotient in C and control is transferred to the divide sub-routine. 
After completion of the division, control is returned to this sub-interpretative routine 
and the last instruction of the divide sub-routine is restored to normal. 

14. Interpretation of 10 

The factor B is obtained and control is transferred to the divide sub-routine. 

15. Interpretation of 11 

The sign of the factor in location A is examined. If the sign is minus, control 
is transferred to location B by modifying the general interpretation routine. Sequencing 
of instructions then continues in the normal fashion beginning with B. If the factor in 
location A is positive or its mantissa is zero, sequencing of instructions continues in 
the normal fashion. 

16. Interpretation of 12 

The general interpretation routine is modified so that the next instruction is 
taken from location A. Sequencing of instructions then continues in the normal fashion 
beginning with A. 

17. Interpretation of 13 

The factor in A is examined and if its mantissa is not zero, the general 
interpretation routine is modified to take the next instruction from location B. Sequencing 
of instructions then continues in the normal fashion beginning with B. If the factor in A 
is zero, sequencing of instructions continues in the normal manner. 

18. Interpretation of 14 

The factor B is obtained, the last instruction in the multiply sub-routine is 
modified to store the product in 0037, which is A, and control is transferred to the 
multiply sub-routine. After completion of the multiplication, control is returned to the 
sub- interpretative routine, the last instruction of the multiply sub-routine is restored 
to normal and control is transferred to the addition sub-routine which computes the 
desired sum. 
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19. Interpretation of 15 

The sign of the factor B is reversed and c ontrol is transferred to the sub- 
interpretative routine "Interpretation of 14". 

20. Interpretation of 16 

This sub-interpretative routine computes the required square root and stores 
it in location B. The initial approximation is X^= (l+4A)/(4+A), obtained from the 
RAND Corporation "Approximations in Numerical Analysis" form 15S, Notes. 

21. Interpretation of 17 

Factor B is obtained and temporarily stored in 0265. Factor A is squared 
and stored in 0057 by the multiplication sub-routine. Factor B is then squared and 
added to the square of A. Factor C is then obtained, squared and added to A* + B*. 
Finally, control is transferred to the square root sub-routine and VA'+B^+C^ 
is computed and stored in 0057 which is K. 

22. Interpretation of 18 

Operation 18 is a vector by vector multiplication. Address Ai and address 
Bx are the addresses of the first elements of each of the vectors. Succeeding elements 
of the vectors are then taken from consecutive memory locations starting with the 
initial locations Aj and Bi- The third address, n, indicates the number of elements in 
each vector. 

Use is made of the sub-interpretative routine "interpretation of 14" to 
perform accumulative multiplication. As each pair of factors are multiplied together 
and added to the previous sum, n is reduced by 1 and zero tested. When n has been 
reduced to zero, the last two factors have been multiplied and the operation is complete. 

Storage _ 

This system requires 46^ storage locations. It uses locations 0000 through 0465. 
Every location within this block is used, thus making it easy to incorporate this program 
with other programs. The translating routine, TRl may be used to translate this 
program to any desired block of memory locations. It should be translated by an 
even amount however to preserve the even-odd conditions. 

■,.■,-,•■,,. pr, fj ■.-^■:::-y/ .-■■> ' ' ■■■■.■'■ ■■( • ■■•■ •■; 

'■■■•-- ;;/■■ .•. '/ - ^ ./-,.... 



■■:: / 
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REMARKS 



Tnp.rpasp i tn i-t-l 



Pir.k up instnip.tinn i+1 

Is. Instruction ]+^ t.n he 
interpreted? 



Put data at address A in location 



nna?. 



Transfer to proper sub-routine. 



Pnngtnnts 
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PROBLEIV 


1: (A) + 


(K)-H 
aiAi 


► K Subroutine 
,:(K)=a2A2, a=x3 


WRITTEN BY: 




c, A= X. XXXXXXX+ 


LOCATION 
INSTRUCTION 


OPERATION 


ADDRESS 


REMARKS 


AB8RV. 


CODE 


DATA 


NSTRUCTION 


0054 


RAL 


65 


0057 


0061 


\ 0057: K 


0061 


SLT 


35 


0002 


0067 


) 


0067 


STTI 


21 


0073 


0027 


1 aj— ♦0073 


0027 


RTT. 


9.0 


0081 


0034 


) A^-*0081 


0034 


RAL 


65 


0037 


0041 


orta?: A 


0041 


ST.T 


35 


0002 


0048 




0048 


STU 


21 


0052 


0055 


a2— ♦0052 


0055 


STT. 


20 


0059 


0062 , 


Ao -*0059 


0062 


RAABL 


67 


8003 


0069 


" 


0069 


SABL 


18 


0073 


0077 


1 ao 1 - 1 a^ 1 . Modify SRD instruction 


0077 


SLT 


35 


0004 


0087 


for case Vhere 0<[a2 - ajl < 10 


0087 


T.n 


fifl 


0090 


0043 




0043 


ST DA 


22 


0098 


0051 




0051 


BRNZ 


45 


0104 


0056 


Is a2 = aj^ ? 


0056 


RAT, 


6S 


0059 


0078 


ao = a^. Thus Ai = Ao + Ai 


0078 


AT, 


15 


0081 


0038 




0038 


BEJaZIL. 


44 


0091 


0044 


Tr lAi i >10? 


0091 


SRD 


31 


0003 


0101 


I |A^|>10 Thus A3 = 0.1 Ai 


0101 


SLT 


35 


0002 


0058 




0058 


RR MIN 


4fi 


0111 


0112 


Tr A3 <0? 


0111 


SIi 


16 


0066 


0070 


Ag-^Oi Tnfi'ea.se exponent hy 1 


0070 


SABL 


18 


0073 


0127 




0127 


SLT 


35 


0008 


0045 


Put exponent in normal positions 


0045 


AT. 


15 


8003 


0053 


(2 hig^h order riigiit positions) 


0112 


AL 


15 


0065 


0120 


A3 •n"'0. Increase exponent by 1 


0120 


AABL 


17 


0073 


0127 




0044 


RRN7. 


45 


0049 


0053 


|A|jK 10. Tr Ai= 0? 


0049 


RAU 


60 


8002 


0066 




0066 


RCT 


3fi 


0000 


0088 


0<fAil<10. Thus A3 equals Ai 


0088 


SU 


11 


8002 


0105 


shifted and counted. 


0105 


RAT, 


65 


8003 


0116 




nilfi 


pp? MTN 


4R 


0070 


0120 


Is A3 ^"0 ? 


0104 


T.Tl 


69 


0057 


0060 




0060 


SRT 


30 


0005 


0074 


• ao/ ai. Is Ia2 - a. 1 >10? 


0074 


RRNZ 


45 


0080 


0079 




0080 


BBJMm_ 


4fi 


0033 


0084 


1 a-j - aj i >10. Is ao >vl-^ ? 


onas 


RAT. 


65 


8001 


0053 


a.*Sao4.10. Thus £3 A3 = ai Aj 


0084 


RAL 


65 


0037 


0053 


ai '^a!t + 10. Thus a^ A3 = a? A« 


007fl 


PR MTN 


46 


0032 


0039 


0V| a\. - fli 1 <r 10." TR"ao ->ai ? 


0032 


RAT. 


65 


0059 


0064 


"\ ' - 1 ' - 1 


0064 


SRT 




0002 


0072 


? ai ">ao^ Thiip Ak = Ai 4- shifted Ao 


0072 


T-P 


69 


0081 


0098 


I X £t *• ** 


0Q98 


SRD 


31 


(0000) 


0071 




0071 


ST.T 


35 


0002 


0030 


Pnmpiite k\ 


0030 


AL 


15 


8001 


0038 




0039 


LD 


69 


0052 


0068 




0068 


STD 


24 


0073 


0076 




0076 


RAL 


65 


0081 


003fi 


an >ai. Thus A| = Ao + shifted AjL 


003R 


SRT 


30 


0002 


0046 


AKsn^'iiRe a2 to coTnpiiVe ag 


0046 


T.n 


69 


0059 


0098 
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PRORI FM; (A) 4- (K) »K Subroutine 



WRITTEN BY: 



LOCATION 

OF 

INSTRUCTION 



00R3 



0090 



OOBS 



(A)=xa. Ai, 



OPE RATI 



TON 



ABBRV. 



STL 



CODE 



^fl- 



31 QQII0Q07L 



nnnniionooL 



(K) a ag A9 a = XX. A a X. xx^ss^xj- 

ADORESS 1 "°" 



DATA instruction! 



nos7 



REMARKS 



00^6 



Store sum in K 



Constants 



PRORl FM. (A) X (B)-*k SubrnutinR WRITTEN BY:_ 

(A) = mi Ml , (B) = m2 M2 , m=xx, M=x.xxxxxxx + 



location 

of 
instruction 



0134 



0141 



0097 



Q1Q6 



0063 



Q119 



J112a_ 



00«R 



0093 



0096 



0103 



0161 



0117 



0075 



0138 



0147 



Q1S5 



0113 



-0121- 



niiR 



0095 



ona?. 



0135 



0086 



0094 



0100 



0107 



01 R3 



0170 



0123 



0181 



-0102. 



OlOfi 



0114 



0171 



0124 



_Q12S_ 



QQ83 



ni9fi 



_J1Q42. 



-JQOOO- 



0050 



0131 



-0139- 



ni7» 



0109 



operation 



ABBRV. 



RALi 



ST.T 



STL- 



BAR. 



SET- 



BAABLl 



SL- 



AIL 



BE-MIN- 



SIL 



SLu 



SLT 



AIL 



RSU 



SBJL 



SIL 



AL_ 



ST.T 



AJJ^ 



RATI 



MBY- 



RR MTN 



SL- 



AL 



IiD_ 



SJLT- 



code 



65 



as. 



20- 



60 



30- 



01. 



16- 



IflL 



40- 



IL 



16- 



25u 



» 



30 



11 



15_ 



^3- 



10- 



fiO- 



10- 



M. 



lfi_ 



15- 



fifl. 



SRIL 

SCT 



RROV 



SRT 



SRIL 



ST DA 



ST TA 



RAI.I 



BR MTN 



AL_ 



SL- 






RAT, 



RRMTN 



Six 



AL- 



-5000001000 



^400000)000 



nnnnnnmoR 



010000)000 



36- 



43- 



M. 



31. 



22- 



23- 



05- 



46- 



15 



IfiL 



20 



65 



4&. 



16- 



15- 



ADORESS 



DATA 



instruction 



0037 



nnna 



J)151- 



.2003- 



0002 



8002 



0131 



0089 



0096 



8001 



8001 



J1D02. 



8002 



8003 



0002 



8001 



8001 



jmQ2- 



-8002- 



8003 



0151 



0135 



0139 



MM. 



8003 



-00 



-00 



& 



0001 



_Q123_ 



QQQ9 



0002 



0205 



0211 



-8001. 



0124 



Q178 



0178 



0Q37 



«oni 



QOOQ 



0109 



0109 



0141 



0097 



0106 



0063 



0119 



0128 



0085 



0093 



remarks 



0037; A- Sfiparatft exponent and 



mantisga 



_Mi=Jt0151 Ml = x . xxxxxxxxoo 



-Jm 



-50- 



0147 



0103 



0161 



0117 



0075 



0138 



0095 



0155 



0113 



0121 



-0115- 



0138 



-0082- 



0086 



0094 



0Q94 



0100 



-0103- 



0163 



0170 



0126 



0181 



0102 



0108 



0114 



0171 



0125 



QQ83 



0083 



QQ26 



0042 



QQ50 



0083 



-0083 — I 



0089; B. Test to srp if M2<Q 



-M2-<J[1L 



Put Mo In iipppr a p.r.inrmlatnr and 



Mq = mj + mo-SO in lower accumulator. 
Adc; odMMMMMMMM/mmoooooooo 



M» > 0. 



Piit M2 in uppftr anp.iimula tnr and 



m^ = mi + m2 - 5 in lower ac c umulatnr 



M§ = Ml X M2 



Round M.^ in 7th decimal place 



Tpwt tr» gpp if M^ is nf thp fnrm 



ox.xxxxxxxx or xx; xxxxxxxs- 



M^ is nf the form ox.xxxxxxxx 
T^niinri M^ In »th rippimal plarp 



to gftt Mg 



'y Mk is nf the fnrm xx. xxxxxx. Thus, 
> mg = mj +1 



t 



Const a nts 



;v- *4m. ^£nJSiy *U-C**m^ „4fUtJ^ , O/Tt .A*rCl/ y'U.^^,iJl^ ^ 



*tr-l^i 






/tf-^u 
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PRDRI FM- (A)/(B) ^^-»K Subi-outine 



(A 



f^ 



diDl 



(B) = dg D2 d=xx, 



WRITTEN BY:. 
D= X.: 



LOCATION 

OF 

INSTRUCTION 



OPERATION 



ABBRV. 



CODE 



ADDRESS 



DATA 



NSTRUCTION 



REMARKS 



0213 



StX- 



55- 



-0002- 



-0169- 



nififl 



RTTT 



9.1 



01 74 



0177 



d2^'»ni74 



0177 



STT. 



20- 



0081 



0184 



D2--»0081. 



0184 



RAL 



fiS 



0037 



0191 



191 



SLT 



5L. 



0002 



0148 



0148 



STU 



21 



00S9 



_Q2a3- 



d^~-»OQ59 



0263 



RAU 



60 



8002 



0221 



0221 



SRT 



5Q_ 



0001 



_Q22!L 



n j) = ni/D2 - 



0227 



DIVR U 



M. 



081 



0464 



0464 



£RJVIIIi 



46_ 



0159 



0209 



0159 



SL 



lfi_ 



0162 



0167 



Round D^ in 8th decimal place. 



02 09 



AL 



Ji 



0162 



_Q1M_ 



0167 



ST.T 



.35- 



0001 



0223 



Test to see if rounded .Do is nf the 



■ J- in III ^.iit 



_Q22a_ 



BRN7TT 



44 



0130 



n9.9.« 



ffirm TT- yyxyyyyyx nr 



-022a_ 



BRMTN 



46- 



0132 



0182 



0132 



AL_ 



15l 



0185 



0189 



n^ is nf the form o. xxxxxxxxx, thus 



0182 



5JL- 



16- 



0185 



0189 



r o v u id Do in the 9th decimal place 



01 89 



SRT 



M- 



0002 



0145 



tn get Dg. 



0145 



SL T 



15_ 



0002 



0152 



0152 



BRMTN 



4fi. 



02 55 



0156 



0255 



SABT. 



18- 



0059 



0164 



0164 



Six. 



1£_ 



0217 



0321 



0321 



AABL 



12_ 



0174 



0179 



d 3 = di - do 4- 50-1 



0156 



A ART, 



12- 



0059 



021 4 



0214 



AIl. 



15- 



0217 



0371 



0371 



MBL_ 



1S_ 
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0241 



0344 



0197 



nns3 



0350- 



3205- 



I 



Restore last instruction of add 



rniiHnp to norma 1 

Prftpare tn compute -^ {fiiYMf^V+i^ - 
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FLOATING DECIMAL SUB-ROUTINES FOR THE IBM TYPE 650 
G. R. Trimble, Jr. 



' Introduction 

The sub-routines given in this article will perform the basic operations of 
addition, multiplication, division and square root using a floating decimal point 
number system. These sub-routines are designed to obtain maximum speed at the 
expense of doubling the memory needed for storage of data. The coder places the 
factors in specified registers, places the next instruction in the distributor and 
transfers control to the sub-routine. The floating decimal operation called for is 
performed, the result is left in the lower accumulator and in the distributor an& 
control is returned to the main routine. Usually, the next operation performed will 
be to store the result in the desired locations. 

The entire set of. sub-routines uses locations 0000 through 0193. It is 

4 programmed to use every memory location within this bloclc of memory and does 

' not use any memory locations outside of it. 

y Each number is composed of 2 parts. The mantissa is 10 digits in length and 

may be either positive or negative. The exponent is 6 digits in length and may be 
either positive or negative. There is no special relation between storage locations 
of the mantis sae and exponents. 

Discussion of Sub-Routines 
> The system consists of the following sub-routines. 

Operation Estimated Average Time 

I' Addition 18.4 ms 

Multiplication 17. 5 ms 

Division 17. 5 ms 

Square Root 118.6 ms 

Subtraction can be performed by simply reversing the sign of one of the 
factors before performing an addition. 

The times given above do not include the time required to obtain factors, place 
them in the calling registers, and store the result. This requires approximately an 
additional 30 ms for the add, multiply and divide sub-routines and 15 ms for the square 
root sub-routine. If results are stored judiciously, it will not always be necessary to 
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obtain both factors. For example, if a product was to be used in the addition sub- 
routine, it could be stored in the calling registers for the addition sub-routine and thus 
already be in place at the time it was needed. Use of techniques such as this can reduce 
the time required for calling sequence by a factor of about 1/2. 

It is assumed that factors are in "normal "form when a sub-routine is entered. 
That is, the high order digit is non-zero unless the entire number is zero. If this is 
the case, the result will always be in normal form. 
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0167- 



.X)122. 



JM.36_ 
.0093 



-0151- 



_DBLE1_M 



..AIl. 



_SEIL, 



.AIL 



.AL 



.AL- 



_AL_ 



RAL 



AU 



.SL_ 



DIVRlL-61 



AL 



SRT 



JLIL. 



AT. 



AL 



_AL- 



RAL 



_A1L 



_aL_ 



DIVRIf 64 



JdjL. 



SRT 



-AIL 



_L5_ 



_21L 



JUL 



JL5_ 



JL5- 



JLi. 



65 



JUL 



16 



15 



_3Q_ 



JJL 



J15_ 



15 



JL5- 



M- 



JLQ- 



_Lfi. 



Ji. 



_3Q- 



JiQfi2u 




First approximation 



X. =( l-h4A) /( 4 +A ) 



0120- 



0178 



0161 
0181 



_Q13a_ 



Ts A ^ = 0? 



-8001- 



-0115- 



_aQOl_LJllM 



0001 



-0185 



0182_ 



-8002- 



8QQ1_ 



-BQQ2. 



JQQ3. 



8002 



0176 



8002 



8001 



8001 



QQOl 



8002 



8001 



8002 



8003 



8002 



0176 



_aQQ2_ 



8001 



8001 



onm 



in I 8002 



J1192_ 



_Q14.9_ 



J115&. 



_Dl-65- 



0123 



0082- 



0092 



-0186- 
0190 



01M_-Tx^ = 1/2 (X^ + A/X^) 



JllM- 



0138 



0147 



0112 



0172 



0132 



0142 



0188 



0140 



-0159 



0169 



X^ = l/2 (X^ + A/X^) 



1 X3 = 1/2 (Xg + A/Xg) 
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IBM 



IBM TYPE 650 PROGRAM SHEET 



FORM NO. 22-6181-0 
PRINTED IN U.S.A. 



PROBI FMr FLOATING DECIMAL SQUARE R OO\R/RITTFM BY:. 



LOCATION 
INSTRUCTION 


OPERATION 


ADDRESS 


REMARKS 


ABBRV. 


CODE 


DATA 


NSTRUCTION 


nifi9 


AT. 


1R 


8001 


0187 




0187 


AL 


15 


8002 


0146 


■ 


0146 


AL 


15 


8003 


0162 


X„ = 1/2 rx, + A/X^) 


8183 


RAL 


65 


8002 


0173 


to £i £i 


0173 


AU 


10 


0176 


0135 


\ 


9M5 


ST, 


16 


800?, 


0144 




0144 


DIVRl 


' e)4 


8001 


0189 




0189 


AL 


15 


8001 


0148 


. 


0148 


SRT 


30 


0001 


0108 


. X/ts 1/2 (Xt 4- A/Xo) 


0108 


AU 


10 


8002 


0168 




0168 


AT. 


15 


8001 


0175 




0175 


AL 


15 


8002 


ooa.'i 




008S 


AT, 


15 


8003 


0143 




0143 


RAU 


fiO 


8002 


01 R9. 


/ 


01 S2 


SCT 


36 


0000 


OlfiO 


f A = V A = X Rhiftpri nr 1 


_aifii_ 


BROV 


47 


0115 


0115 


\ £i X '±. 


0115 


STU 


21 


0170 


0124 


\ Put a, in lower accumulator, A„ in 


0124 


RAL 


65 


0177 


0184 


1 distribuior and go to next inst. " 


0184 


T,n 


69 


0170 


nnfl7 








































0107: 


0100000 


non 






^ — . 


0167: 


0400000 


000 






Constants 


0104: 


0000500 


UOO 
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IBM TYPE 650 LOADING ROUTINES 

G. R. Trimble, Jr. 
E. C. Kuble 



One of the first problems that arises in using a computer of any sort is how to 
enter information into the machine. The flexibility of the input system of the IBM 
Tjrpe 650 gives the programmer a wide range of possibilities. Most problems will 
require their own individual loading routines which are adapted to the particular 
type of data which they will be processing. The following routines, however, are of 
general interest since all problems require that programs and possibly tables, be 
entered. These routines are by no means exhaustive but merely indicate the many 
possibilities which exist and the ease with which this problem may be solved for the 
650. 

These routines are programmed so that the card reader will operate at 200 cards 
per minute. In some cases, it was necessary to use optimum programming to obtain 
this speed. In addition, they are programmed so as to use a block of memory 
location. The translating routine may be used on each of these programs, as well as 
the translating routine itself, to store them in a different set of memory locations from 
those for which they are programmed. 

Since unpunched columns enter the machine as. invalid information, the word entered 
will.notpass the validity check if the program tries to use it. By punching zeros with a 
plus sign in all unused fields of the different types of loading cards, the information which 
enters the machine will be valid though perhaps meaningless. Thus, if it is necessary 
to pimch out the contents of the entire memory in order that a new program can be 
entered, it is not necessary to make provision for avoiding locations in which invalid 
information may possibly have entered. This procedure greatly simplifies the routine 
which would be required for unloading memory. 

DESCRIPTION AND FUNCTION 
LLl Loading Routine 

This routine enters four words from a card, each card containing the addresses of 
the locations in which these four words are to be stored. It consists of nine instructions, 
five of which are kept on the drum in locations 1995 through 1999, the other four being 
entered on each card into locations 1951, 1953, 1955, and 1957. The four words to be 
stored are entered into locations 1952, 1954, 1956, and 1958. The routine consists 
essentially of a read instruction followed by four sets of load distributor- store 
distributor instructions. The data addresses of the store distributor instructions, 
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which are punched on the load card, ' specify the locations into which the corresponding 
words are to be entered. Load cards are identified by a 12 punch in column 1. The 
only control panel wiring necessary is from column 1 of first reading to the load hub. 

LLl is particularly useful for entering instructions. The card form is designed 
so that the card can be punched directly from the program sheet. The address of the 
location to be entered is punched followed by the word itself. 

This routine is also useful for entering programs which are coded optimumly where 
instructions are scattered throughout the drum. It can also be used to load other loading 
routines which are not self-loading. 

LLl can be loaded with only two cards. To use LLl simply place these two cards 
in front of the program to be loaded, set 70 1951 3000 on the storage entry switches and 
take the first instruction from the storage entry switches. This causes cards 1 and 2 
to read into the 650 and the program punched on these cards enter LLl in the desired 
locations. Control is then transferred to LLl and the cards following are entered by it. 

Once all of the desired cards have been read, it is necessary to transfer control 
to the main program. This can be accomplished in one of two ways. If the card read 
is not a load card, the next instruction will not be taken from 1995 as is usually the 
case, but will be taken from the location specified by the I address of the read 
instruction. Thus, the I address may be changed to cause control to transfer to the 
first instruction of the routine when the program has been entered. 

The second and probably most useful method is simply to change the I address 
of one of the store distributor instructions punched on the load card. Thus, if only 
two instructions are punched on a particular load card, the I address of the second 
store distributor instruction, which is in columns 27 through 30 should be changed 
to cause control to transfer to the desired location. Obviously load cards containing 
one, three, or four words can be handled in the same manner. 

In case only one, two, or three words are entered from the card, the remaining 
fields should have zeros punched. This makes it easy to punch out the contents of the 
entire memory at a later time if desired. 
Load Card Form - 12 punches in columns 1, 10, 30, 50, 70. Signs of words 1, 2, 3, 



and 4 must be either 11 or 12 punch in columns 20, 40, 60, and 80 respectively 



IfSX 



.1 

n 

0000 

3 4 3 S 

nil 

2222 

3333 
4444 
5555 
S66E 
7777 
888S 
9999 



0000 

7 8 9 10 
gill 

2222 
3333 
4444 
5555 
656E 
7777 



9|fl9 
T e 9 i( 
5oal 



iq&z 



Word 

0000000000 

r 12 13 14 15 IS 17 IS 19 20 
1111111111 

2222222222 
3333333333 

4444444444 
5555555555 
G 630666656 
7 7 77 7 77777 
88888686 88 
9999999993 

11 13 13 14 IS IS 17 IS 19 20 



iqffa 



.5 wj 
5 3 

0000 

23 24 35 26 
1111 

2222 
3333 

4444 
5555 
686E 
7777 
888S 
9999 

;4 2S2E 



000( 

27 20 19 3 

Bt 1 1 

222! 
3333 
4444 
5555 
6666 
7771 
8888 



9IB 

27 29 30: 



>qffM 



Word a 

0000000000 

31 33 33 34 3S 30 37 30 
1111111111 

2222222222 
3333333333 
44 44444444 
5 5 5 5 5 5 5 5 5 5 
6606666668 
7777777777 
8088888888 



»9 



999999993 

32 33 34 35 30 37 30 39 



(isy 



S -o 



0000 

43 44 45 40 
1111 

2222 
3333 
4444 
5555 
6666 
7777 
8888 
9999 

43 44 45 40 



I9ff6 



Word 3 



0000000000 

pi 52 S3 54 5£ so 57 50 59 OO 
1111111111 

2222222222 
3333333333 
4444444444 
S5SSS5555S 
6G66G6G86E 
7777777777 
388888888! 
9993999939 

47 % 49 S( >l 52 53 54 55 56 57 50 SO CI 



0000 

17 40 49 50 
1111 

2222 
3333 
4444 
5555 
6666 
7777 
888t 

f> 

140 491 



fTsr 



c 

0000 

13 U O'l EG 
1111 

2222 
3333 
4444 
5555 
6666 
7777 
8888 
9 9 9 9 

El 62 03 04 65 66 



OGO 

67 60 63 

B11 

222 
333 
444 
555 
666 
777 
888 
98(1 

67 60 69 



msa 



Word f 

0000000000 

71 72 73 74 75 76 77 73 79 00 
1111111111 

2222222222 
33333 33333 
4444444444 
5555555555 
6666666666 
7777777777 
8888888888 



9999999999 

7d n 73 73 74 75 70 77 70 19 SO 
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IBM 



IBM TYPE 650 PROGRAM SHEET 



FORM NO. 22-6181-0 
PRINTED IN U.S.A. 



PRHRI FM. LLl LOADING ROUTINE 



WRITTEN BY:. 



LOCATION 
INSTRUCTION 



1995 



1B2£. 



1997 



199R 



1999 



TMT 



1953 



1955 



1957 



OPERATION 



ABBRV. 



LP 



ML 



Ul_ 



UDl 



HEAD. 



STD 



ST D 



STD 



STD 



PKpm.T IM- TnADTNf BOTTTTT^K TO T.f AD-LLl 



Bonn 



1951 



1901 



1902 



1903 



lSft4_ 



1905 



iaofi_ 



1907 



1957 



1958 



HEATL 



READ 



LIL 



STD 



CODE 



_fiS_ 



m. 



-E9- 



JZO- 



24 



J24_ 



_24. 



_24_ 



-70- 



JUL 



_&9. 



LD 



STD 



LEL 



STD 



LIL 



liD- 



STn 



24 



ADDRESS 



DATA INSTRUCTION 



-ISSi 



1954 



195fi 



1958 



1995 



xxxx 



xxxx 



xxxx 



xxxx 



1951 



_19Qi 



1951 ^ 



1953 



1955 



-1.95!3^ 



0000 



1996 



1997 



199R 



1999 



3000 



1952 



69 



_24. 



24 



69 



_24. 



-sa. 



.^4_ 



1995. 



1953 



3noo 



1902 



1903 



IQQfi 



1954 



1997 



1955 



1998 



1956 



1999 



1 904 



19n5 



REMARKS 



These instructions are kept on the drum. 



Thp.dP ingtrnptinn jg arp r pad into the 650 



on each load card. The data addresses 



specif y the Inp.atinn where each of the 



) four pieces of information punched on 
the card are to be stored. 



Read Card 1 intn 1951 thru 1958 



\ Read Card 2 into l a Ol thru 1908 



1906 



19Q7 



1908 



1957 



1958 



19Q9 



I 



f3f,^no T.T.I in iggR thru 1999 



1952: 



1953r 



1954; 
1956; 



69 1952 



-195; 



69 1954 



jLaai 



69 1956 

69 1958 

70 1995 



195Si 
000(p 



PAGE. 



.OF. 



46 



Cards 1 and 2 for LLl 

Card 1 is punched as follows: 12 punches in columns 1, 10, 20, 30, 40, 50, 60, 70, 80 

101901300Q g91952195ia919541953 &919 561955691958195Z701995000Qg9195619512419991 99§, 
1951 1952 1953 1954 1955 1956 1957 1958 

Card 2 is punched as follows: 12 punches in columns 1, 10, 20, 30, 40, 50, 60, 70, 80 



^9i95 2l90^g419 95190j^g919531904 2 41996190^g9195419062419971902 £91955190§ g41998195X 
1901 I9S2 1903 1904 1905 1906 1907 1908 



LI Loading Routine 

LI is used to enter load cards as identified by a 12 punch in column 1. As many 
as seven ten-digit words may be entered on a card by means of LI. A control word 
punched on the card indicates where the first word is to be entered as well as the 
number of words to be entered. Successive words on the card are entered into 
successive memory locations. 

Card Form for LI Load Card: 



1 10 


11 20 


21 30 


") 


Control 
Word 


Wordl 


Word 2 


Wo/ 



70 71 



80 



Word 7 



The control word is constructed as follows: 

Column 1, double punch 12 and 0; 

Column 1, punch 0; 

Columns 3-6, punch address where Word 1 is to be entered. 

This address is f=f 1 f^ isU. 

Columns 7-9, punch 0; 

Column 10, double punch 12 and n, where n is number of words to be entered. (n=l, 

2, 3, 4, 5, 6 or 7) 

Load cards must have all ten columns of the fields to be entered punched and must 
have the sign punched over the units position of the field (columns 10, 20, 30, etc.) 
12 for plus, and 11 for minus. Fields not entered should have zeros punched in them. 
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Control Panel Wiring: 

Column 1 of first reading is wired to the load hub. The 12 in column 1 identifies 
load cards. 

Explanation of Program: 

Instruction 1988 calls for a Read. K the card read is not a load card, the data is 
entered into locations 1951-1960 and the next instruction is taken from m, where m 
is the location in the main routine to which control is transferred. If the card read 
is a load card the data read is entered into locations 1951-1958 and the next instruction 
is taken from 1994. 

Words 1 through n will be transferred to locations f through f + n - 1 by LD and 
ST D instructions. These two instructions will be kept in the accumulator where they 
will be modified and executed. 

Instruction 1994 puts the control word in the accumulator. 

Instructions 1979 and 1986 modify the ST D instruction so that Word 1 is stored 
in location f. The ST D instruction is then 24 (f) 1977. 

Instructions 1989, 1981, and 1990 construct a dummy instruction which is used 
to indicate that the last word has been transferred. This dummy instruction is 
24 (f + n) 1977. 

Instructions 1984 and 1992 load the accumulator with the LD and ST D instructions. 
These instructions are then executed from the accumulator. 

Instructions 1977 and 1985 add 1 to the data addresses of the LD and ST D in- 
structions in the accumulator so that the next time they are executed, the next word 
will be transferred to its proper location. 

Instructions 1993 and 1983 compare the ST D instruction in the accumulator with 
the dummy ST D instruction. If they are equal, as indicated by a zero in the upper 
half of the accumulator, the last word has been transferred and control is returned 
to 1988 and a new card is read. If they are not equal, the ST D instruction is 
regenerated in the upper half of the accumulator by adding the dummy instruction 
back in on instruction 1987. Control is returned to 8002 and the next word is 
transferred. 
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IBM IBM TYPE 650 PROGRAM SHEET 

irwilcMntV 

PRORI FM. LI LOADING ROUT INE WRITTEN BY: 



FORM KO. 2Z-6ISI-0 
PniNTED IN U.S.H. 



LOCATION 

OF 

INSTRUCTION 


OPERATION 


ADDRESS 




aeu&auti 


ABBRV. 


COOE 


DATA 


NSTRUCTION 


'"• 1 


1^88 


READ 


70 


1994 


m 




(Transfer control to location m If no load card) 


1994 


RAL 


AS 


1951 


1979 






1979 


LD 


69 


1982 


1986 






19R6 


ST DA 


22 


1982 


19ft9 






1989 


SLT 


35 


0004 


1981 




^Propnm tn ^ynra «tnrtlng with location f . 


1981 


AL 


15 


8001 


1990 






1990 


ST DA 


22 


1978 


1984 






1984 


RAL 


65 


1991 


1992 






1992 


AU 


10 


1982 


__JB0D2 






Rnn? 


in 


69 


(1952) 


8003 


\ Store word [ in location f + | - 1 


8003 


ST D 


24 


(f) 


1977 




1977 


All 


10 


1980 


1985 


\^AJdIfy Instructions In preparation to store 


1985 


AL 


15 


finoi 


1993 


/ viorA '[ + 1 . 


1993 


SU 


11 


1978 


1983 


"\Test to see If word n - 1 has been stored. If 


1983 


RRNZtl 


44^ 


1987 


1988 


/so, read the next card. 


1987 


AU 


10 


8001 


_8002- 


y Propnm in stnrfi word j + 1 - 


f 












1982 


STD 


24 


QDQil- 


1977 


] 


t99T 


LD 


69 


1952 


8003 


^Constonts 


T980 




00 


0001 


0000 


1 














1978 


STD 


24 


f +n 


1977 


^ Temporary Storagfi 














1951 




00 


f 


OOOn 


y Control word punched on load card . f Is 












flddreSS whnrfi f rst word 0" r.ord s tn hn 












entered . n Is t ie number of words to be 












entered. 
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L2, L2A, L2B, L2C. Loading Routine 

L2 is used to load eight full words from each card. The initial address f^ is 
given, and words are located sequentially beginning at f . Cards containing a 12 in 
column 1 will be loaded in consecutive addresses. This will continue until occurrence 
of a card with no 12 in the load control column at which time automatic branching is 
made from L2. 

An optional way to get out of the routine is to set the address switches at the first 
location beyond the last one for which loading is desired. Then branching can be 
manually effected to other routines. In this case the word following the last word 
to be loaded must have ten zeros punched. 

The calling for L2 and f can be effected manually, from the first card read, or 
from the preceding program by L2A, L2B and L2C respectively. 

Wiring: 

Column 1 of first reading is wired to the load hub. The 12 punch in column 1 
identifies load cards. 

Program Explanation: 

The store distributor operations are executed from the lower accumulator and 
stepped up by the contents of the upper accumulator. The load distributor operations 
are stored in memory. Loading is then realized by repeated load distributor, store 
distributor (per instruction in lower accumulator), and add to lower accumulator from 
upper operations. 

L2A: L2A is designed to initiate L2 and specify f^ manually. The storage entry 
switches are set to 00 f 1966 and control is sent to these switches by depressing in 
turn the program reset button and the program start button. L2A then inserts f^ 
into the L2 routine and reading is begun. 

L2B: L2B is designed to have the first card of a group to be loaded by L2 designate 

f . The first word on this card is punched 24 f„ 1981. L2B transfers f„ into L2 and 
o o o 

initiates L2. To call for L2B all that is required is to give an instruction address of 
1970. 

L2C: L2C is designed to have f designated by a word in memory. This word is in the 
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following form; 24 i^ 1979. L2C inserts f^ into L2 and initiates L2, To call for L2C all 
that is required is to give an instruction of the form 65 m' 1969. 

Fields on the last card which are not entered should have zeros punched in them. 



IBiyi IBM TYPE 650 PROGRAM SHEET 

Tr;iac-Mark 

PPmi FM. T.9.A^ P>, and C. T-OADTNG ROUT INEWRITTEN BY: 



FORM NO. 22-6I(!t-0 
PRINTED IN US A. 



LOCATION 

OF 

INSTRUCTION 


OPERATION J 


ADDRESS 


REMARKS 


ABBRV. 


CODE 


DATA NSTRUCTION 










i 


L.2 A (fr 


takpn f 


•nm p!' 


nrntrp Krl*"'''^ Sij7itfihpR\ ... 












8000 


No Od 


GO 


f 


1961 


14fiR 


T?AT, 


65 


fioon 


1 974 
__ia6.2.-.i 


1 


1974 


SRT 


30 


0004 


_ ( 


1967 


SLT 


35 


0004 


1968 i 
196a^J 




1968 


AL 


15 


1971 "^ 


L. Ace: 24 fo 1979 


19R9 


ATT 


10 


1972 


„ 19.11- ; 


U. . Acci Qft-O'QQl 0002. _ i 










i 


1P71 


pT n 


24 


onnn 


1979-^ 


[ Constants ... . . i 


1 Q79 




Od 


nnni 


0002 


f i 
























■ 












1 














T,5).R (f„ 


tijkfin frc 


km fir 


=;t word o 


f first ca, 


rd) ; 














1970 


READ 


-70 


1975 


0000 




1975 


RAL 


65 


1951 


1976 


L. Ace: 24 f^ 1981 


1976 


ATI 


10 


1972 


1980 


U. Ace: 00 O^iOl 0002 ! 












-.---. ^ . , - - i 


IPFil 


RT r> 


24 


ffi 


1981 


T?irst word punched on first card | 














1972 




no 


0001 


0002 


Constant 




































1 














T.2C: (fr. 


in memo 


rv loc 


ation m') 


















m" 


RAL 


65 


m' 


1969 


L. Ace: 24 f_ 1979 


IQRS 


ATT 


10 


1972 


i 1011 


TT. Ace: 00 0X)01 0002 










i ' 


m' 


ST D 




24 


^o 


1979 i 




















































































1 




















































































! 












1 
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IBIVI IBM TYPE 650 PROGRAM SHEET 

Tmle-Marli 

PRORI FM- L2 LOADING ROUTINE WRITTEN BY: 



FORM NO. 2Z-6ial-0 
PRrNTED IN U.3.A. 



LOCATION 
IN3T8uGTIO^ 


OPERATION 


ADDRESS 




aeutkoieti 


ABBRV. 


CODE 


DATA 


NSTRUCTION 




1977 


RFAH 


70 


1978 


m 


(Jmnttar rnnfrni fn \neat\nn rtl If no load card) 


1978 


LD 


69 


1951 


8002 






8002 


STD 


24 


(n 


1979 






11979 


Al 


15 


8003 


1980 






1980 


LD 


69 


1952 


8002 






8002 


STD 


74 


ff + n 


1981 






1981 


AL 


15 


8003 


1982 






1982 


LD 


69 


1953 


finn9 






8002 


STD 


24 


(f + 9) 


1983 






1983 


Al 


15 


8003 


1984 






1984 


LD 


69 


1954 


8002 






8002 


STD 


24 


(f + 3) 


1985 






1985 


AL 


15 


8003 


1986 




Y Store 8 Wrtrds In locations f thru f + 7. 


1986 


LD 


69 


1955 


8007 






8002 


STD 


24 


(f + 4) 


-J287. 






19S7 


AI 


1,5 


8003 


1988 






T98» 


LD 


69 


1956 


8002 






sma 


ST D 


24 


ft + 5) 


1989 






19^ 


AL 


15 


8003 


1990 






1990 


LD 


69 


19V 


8009 






8002 


STD 


94 


— XLt-6L 


1991 






1??I 


AL 


15 


8003 


1992 






{9^2 


LD 


A9 


1258 


8002 






8002 


STD 


24 


(ft 7) 


1993 






1993 


AI 


15 


8003 


1994 


"^ Modify inslruf^tlon"! to rppnnt Innp for 


1994 


SL 


16 


1973 


1977 


r next card . 














1973 




00 


0000 


0016 


Constant 
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LT 1 Loading Routine For Tables 

LT 1 is used to enter table cards as identified by a 12 punch in column 1. One 
argument and as many as six functions associated with that argument may be entered 
on a card by means of LT 1. A control word punched on the card indicates where the 
argument is to be entered, the number of arguments and functions to be entered, and 
where the functions will be placed relative to the argument. 



Card Form for LT 1 Table Loading Card: 



10 

Control 
Word 



11 



20 



Argument 



21 



30 



First 
Function 




The control word is constructed as follows: 

Column 1, double punch 12 and 0. 

Column 2, punch n, where n is the number of arguments and functions to be entered. 

For one argument n=l, for one argument and one function n=2, etc. n=7 is maximum. 

Columns 3-6, punch address where argument is to be entered. This address is a=ai aa 

^3 a* . 

Columns 7-10, punch increment between argument and first function. This increment 

is i, i=ii is is i4 . Thus the argument will be stored in a, the first fimction stored in 

a + i, the second function stored in a + 2i, etc. Column 10 must also have a 12 punch. 

Table load cards must have all ten columns of the fields to be entered punched and 
must have the sign punched over the imits position of the field (columns 10, 20, 30, etc.) 
12 for plus and 11 for minus. Fields not entered should have zeros punched in them. 

Control Panel Wiring: 

Column 1 of first reading is wired to the load hub. The 12 in column 1 identifies 
table load cards. 

Explanation of Program: 

Instruction 1978 calls for a Read. If the card read is not a table load card, the 
data is entered into locations 1951-1960 and the next instruction is taken from m, 
where m is the location in the main routine to which control is transferred. If the 
card read is a table load card the data read is entered into locations 1951-1958 and 
the next instruction is taken from 1994. 
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Instruction 1994 places the control word in the lower accumulator. 

Instructions 1985 and 1988 modify the ST D instruction so that the argument is stored 
in location a. 

Instructions 1989, 1981, and 1991 store the increment i in 1975. 

Instructions 1984, 1992, 1976 and 1986 construct a dummy instruction which is 
used to indicate that the last word has been transferred. This dummy instruction is 
69 (1953-Hi) 8002. When the LD instruction has been modified so that it is equal to the 
dummy instruction the last fimction has been transferred. 

Instructions 1993 and 1980 load the accumulator with the LD and ST D instructions. 
These instructions are then executed from the accumulator. 

Instructions 1972 and 1979 modify the LD and ST D instructions in the accumulator 
so that the next time they are executed the next piece of data will be transferred to its 
proper location. 

Instructions 1987 and 1973 compare the LD instruction in the accumulator with the 
dummy LD instruction. If they are equal, as indicated by a zero in the upper half of 
the accumulator, the last word has been transferred and control is returned to 1978 and 
a new card is read. If they are not equal, the LD instruction is regenerated in the upper 
half of the accumulator by adding the dummy instruction back in on instruction 1977. 
Control is returned to 8003 and the next word is transferred. 

Note that if i=l, LTl is functionally identical with LI. This increased flexibility 
is gained by an increase in the number of storage locations required however. LTl 
uses 5 more storage locations than LI. 
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IBI^ IBM TYPE 650 PROGRAM SHEET 

Tri.il*. Mirk 

PPnRl FM. LT 1 LOADING ROUTINE FOR TABL ES WRITTEN BY:. 



FORM NO 22-eiai-o 

PRINTED IN USA. 



LOCATION 
INSTRUCTION 



L_i58a_ 

1989 



1^78" 



1994 



JL9-85_ 



_198J_ 



1991 



1584_ 



1992. 



1976 



J3S6- 



J222L 



J58flL 



8003 



8002 



1972 



1979 



1987 



1973 



1977 



19 74 



1983 



1982 



1990 



1951 



OPERATION 



ABBftV. 



READ 



JAL 



JJ3l 



-SJDA 



LD_ 



SLT 



JTDA 



RAL 



SLT 



AL 



_STLDA_ 



RAL. 



JAL 



_LDl 



-STLD- 



AL 



-ALL 



SU 



J&HZU 



ML 



STD 



LD 



JLEL 



CODE 



70 



JiSu 



J&2- 



.21. 



J>9_ 



JJ5. 



J22_ 



.65. 



Ji 



15 



.22. 



jS5. 



_10_ 



J2A. 



Ji 



JjQ. 



11 



U 



JjQ. 



_24 



69 



JIQ. 



ADDRESS 



DATA INSTRUCTION 



"199?" 



1951 



J9^ 
1974 



JI.QQ3_ 



J]D.Q4_ 



J9Z5_ 



_8003_ 

_0002_ 

1983 



J59Q_ 



15Z4_ 



_I98a 

Xi25: 

_L925 
_1582 



12) ^1 



.1?.?D_ 



_15Z7_ 



_aaQL 



.62. 



JMCL 
1952 

JQ-Q-OL 



(1952M 



m 



JS35- 



Jisas. 



1989. 
_198J_ 

_L?84. 
1992 



±976_ 
1986 



1993_ 
.IS8Q_ 
.8003„ 
.8002_ 
_J.9Z2. 
1?Z2_ 



. "\<;trirp wnrfl j in location fl -<• {|-Ili 



_1587_ 

J92a_ 



.L9Z8_ 



_8.aaa_ 



_ISZ2_ 



_8002_ 

_ooao_ 



_8002_ 



REMARKS 

(Transfer control io lo caHon m If no load card) 



i£repgig-tQ-5ioc9 sl grting-Mlt h iocatlon a. 



\ Modify In! 
/ word i + 1 



JnstnjpH""s In prepgroMon to store 



Jest to see if wo 

so, raad the next card. 



Eisjms to store word j + 1 . 



I 



Constants 



J:JEemporflry_Siorage_ 



anitoLwor 

q Is addr ess Where the first word 

■itnrfi d , I Is Incramnnt to be added to addresses 



of siicc»« 1vf! words, n Is the number of words 



JaJjfijstQxesL. 
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TR 1 Translating Routine 

This routine is used to translate a sub-routine or sub-program to a different set of 
memory locations from those for which it was programmed. It will simultaneously 
modify each instruction as required, so that they will be properly executed from their 
new locations. Each instruction or constant in the sub-routine will have associated 
with it a six digit number. Four digits of this number specify its location within the 
sub-routine . The remaining two digits indicate whether or not either the D or I 
address parts of the number should be modified when it is translated. An 8 indicates that 
the corresponding address Is to be modified, and a 9 Indicates that the address Is not to 
be modified. Thus, for example, the modification control digits for a constant would 
be 99 since nothing in the constant is to be modified. Most Instructions, however, will 
have 88 as the modification control digits since both the data and instruction addresses 
will be modified. A shift instruction would usually have 98 as its modification control 
digit, since the data address of a shift instruction should not be modified. Similarly, 
an Instruction in which the D address should be modified while the I address should 
not be modified will have 89 as its modification control digits. 

A master card placed in front of the deck containing the sub-routine indicates the 
amount by which the sub-routine is to be translated, the number of instruction" in that 
sub-routine, and an indication of what should be done after the sub-routine has been 
completely translated. 

The detail cards will each contain five ten-digit words and their corresponding 
control information. If less than five words are punched in the detail card, the 
remaining fields should have zeros punched in them. 

Master Card Form: 

12 punches in columns 1, 10, 20, 30, 40, 50, 60, 70, 80. 

OOOOObttttOOOOOOlnnnnboOOOOliliOO . . .0, 



Amount Number Next 50 zeros 

of of Inst. 

Translation Words 



Detail Card Form: 

The detail card contains five sixteen-digit fields each punched as follo\frs : 
Col. 1: D modification control column: 8 if _data address is to be modified , 9 if data 
address is not to be modified. 



(I 



56 



Col. 2: I modification control column; 8 if instruction address is to be modified, 9 if 

instruction address is not to be modified. 

Col. 3-6: Location of word within sub-routine before translation. Column 6 also has 

a 12 punch. 

Col. 7-16: Punch the. word to be translated. The sign of the word is in column 16 

(units position). 

Control Panel Wiring: 

Wire from column 1 of First Reading to the Load Hub. A 12 punch in column 1 
identifies master cards. 

The six control digits for the first word are wired to enter the 6 low order 
positions of Storage Entry Word 1. Word 1 is also wired for Word Size 6. The first 
word is wired to enter Storage Entry Word 2. The Sign over Units switch and Read 
Plus switch must also be wired. Word 2 must also be wired for Word Size 10. The 
remaining 4 fields on the card are wired similarly to enter Storage Entry Words 3 
through 10. 

Explanation of Program: 

The instruction in location 1933 calls for a read. If the card read is a master 
card, the next instruction will be taken from location 1991 and the nine instructions 
starting with 1991 and ending with 1990 will be executed. These instructions store 
the amount of translation, the number of instructions and the location of the next 
instruction, i, in their proper locations. It places the amount of translation, t, in the 
data address positions of one word, the instruction address positions of another word, 
and in both data and instruction address positions of a third word. Thus, when it has 
been determined how the word to be translated should be modified, all that must be 
done is to add one of these three numbers or zero to that word before storing it in its 
new location. 

The location of the next instruction, i, is used to indicate what should be done 
after the sub-routine has been translated. For example, if another sub-routine foUows 
which must also be translated, i will take the form 00 0000 1933. When the translation 
of the first sub-routine is complete, another read instruction is given to call in the 
master card for the following sub-routine. When the last sub-routine has been 
translated, i is used to transfer control to the first instruction of the program. 

The four instructions, starting with the one in location 1939, simply set two 
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other instructions to their initial values. 

The five instructions beginning with the one in location 1937, add t to the sub- 
routine location of the word being translated and modifies a "store" instruction so 
that this word will be stored in location L + t. Beginning with the instruction in 
location 1968, the modification control digits are analyzed to determine whether the 
data address and/or the instruction address positions of the word to be stored should 
be modified. This is accomplished by adding one of the four constants computed from 
the master card. The instruction in location 1965 then stores the modified word in its 
new location. 

The number of instructions, n, is decreased by one each time a word is stored 
and is zero tested to see if the last word in the sub-routine has been stored. If so, 
control is transferred to i. If the last word has not been stored, a test is made then 
to determine whether the last (5th) word on the card has been stored. If so, a new 
card is read. If more words remain to be transferred, the instructions are modified 
in preparation for storing the next word. 

Use of TR 1: 

Once it has been determined which sub-routines are to be used and where they 
are to be placed, the programmer punches the amount of translation necessary to 
place the library sub-routine in the desired locations and the other information 
necessary, on a master card. This master card is placed in front of the library sub- 
routine and all of the sub-routines are assembled. These sub-routines, along with 
the rest of the program, are then read into the 650 and the problem is begun. 

Routines which are programmed optimumly should be translated only by an 
even amount in order to preserve the even-odd conditions. 
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igj^ IBM TYPE 650 PROGRAM SHEET 

T.i..lr-Mnt>l 

PRORI FM. TR1 TRANSLATING ROUTIN E WRITTEN BY:. 



FORM NO. 22-eiDl-O 
PRINTED INU.B.A. 



. ''"' /^ 



LOCATION 


OPERATION 


ADDRESS 




REMARKS 




INSTRUCTION 


ABBRV. 


CODE 


DATA 


NSTRUCTION 




933 


READ 


70 


1991 


1939 






991 


RAL 


65 


1951 


1938 






938 


ST D 


24 


1992 ^ 


1967 


Sfnm 00 0000 (H 1 




1 


967 


SLT 


35 


0004 


1934 






1 


934 


ST DA 


22 


1941-' 


1929 


Stom 00 (t) (t) 


Master 


1 


929 


ST L 


20 


1993^ 


1940 


Storfl 00 (t) 0000 


\ , Card 


— . 


940 


LO 


69 


1953 


1969 




Only 




969 


STD 


24 


1948-- 


1930 


Store I 






930 


RAL 


65 


1952 


1990 






1990 


STL 


20 


1989- 


1933 


Store n J 




1939 


LD 


69 


1947 


1945 






1945 


ST D 


74 


1949 


1977 




. <;Bt in^triictinn<; 19'^T rinA 1949 tn 


1972 


LD 


69 


1932 


1935 




initlnl vnliiOB. 


1935 


ST D 


24 


1937 


Rnm 


J 




1937 


RAL 


65 


(19511 


1964 


1 




1964 


SLT 


35 


0004 


1976 




Prepare to store a word In location L + t 


1976 


AL 


15 


1993 


1947 




^ 


1947 


in 


69 


19sn 


19A1 






961 


ST DA 


22_ 


1965 


1968 








968 


LD 


69 


ftnn7 


1975 








975 


BRD? 


99 


1979 


1987 




Tflst to sfifi how word should b" 


1979 


PRO 10 


90 


1984 


1988 




mnrlifinH hofi-irn hping stnmd 


19fl7 


RRDIO 


90 


1987 


1986 






1986 


RAL 


65 


1994 - 


' 1949 


1 




1987 


^^L 


65 


1993 


1949 






1984 


RAL 


A5 


1941 


1949 




Modify word and store It In location. L+i.. 


1988 


RAI 


^5 


1997 


1949 






1949 


AL 


15 


(1957) 


19A5 






965 


ST L 


20 


(L + t) 


1974 






974 


RSL 


66 


1981 


1936 


. 


Test to ^pi^ if It* word in rniitinn has he en 


1936 


AL 


15 


1989 


1943 




■ stored . 


1943 


BRNZ 


45 


1946 


1948 




If so, transfer control to I 


1946 


AU 


10 


1949 


1963 






963 


SU 


n 


1966 


1971 




>- Tput tn «PR if Inst word nn rnrd 


971 


BRNZU 


44 


1978 


1990 


) hns hppn stnrRcl. - -' 


978 


AU 


10 


1931 


1985 








J 985 


ST L 


20 


1989 


1944 




■ Prepare to Store next word 




944 


ST 1) 


91 


1949 


19A? 






1 


962 


RAL 


65 


8003 


1970 




' 




970 


SL 


1-^ 


197.^ 


1977 


. Prntwrn tn rantJ nnxtcarJ 




1977 


LD 


69 


1980 


1983 






1983 


ST DA 


22 


1937 


8001 




1 












1 


1941 




00 


t 


t 






J 992 




GO 


0000 

t 


t 








00-5 




00 


0000 




Constnnts determined from MasterCard 




1 


994 




00 


0000 


onon 






1989 




00 


onnn 








1948 




00 


0000 


1 
















1 



1942 




15 


1952 


1965 


"^ 




1932 




65 


1951 


1964 






1950. 




20 


0000 


1974 






1981 




nn 


0000 


0001 




Constants 


19AA 




15 


1960 


1965 






1931 




15 


1962 


1965 






1973 




00 


0001 


0000 






1980 




65 


oonn 


19A4 
















1 
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A METHOD FOR PERFORMING DOUBLE PRECISION ARITHMETIC ON THE IBM TYPE 650 

G. R. Trimble, Jr. 

Introduction 

It is sometimes necessary to perform computations using a larger number size 
than is basically provided on the 650. This system has been developed to perform the 
basic operations of addition, subtraction, multiplication and division using numbers 
which are 20 decimal digits in length. It was designed with coding convenience as a 
prime objective and uses an interpretative routine to perform double precision 
arithmetic. 

Double precision instructions, that is, instructions which are to be interpreted 
have a negative sign. A double precision instruction consists of a 2 digit operation 
code, a 4 digit address specifying the location of the first factor and another 4 digit 
address specifying the location of the second factor. The result of the operation is 
always stored in a pair of memory locations referred to as C. A 5th operation, "Store", 
is used to transfer the result from C to the locations desired by the programmer. In 
order to store the result where desired, it will Xisually be necessary to follow one of the 
4 basic arithmetic operations by a store instruction. The net effect of such a procedure 
is four 3-address operations. 

Tae address given refers to the 10 high order digits of the factor to be operated 
upon. The interpretative routine automatically selects the 10 low order digits from 
the next memory location. Thus, the high order and low order digits of a number will 
always be stored in successive memory locations. 

Should a positive instruction appear, it is interpreted as a normal 650 instruction 
and subsequent instructions are not interpreted. Thus, upon occurrence of a 
positive instruction, the 650 operates in its usual D-I mode. This continues until an 
instruction address of 0004 is given which causes the control to return to the 
interpretative routine. The program will return to the double precision mode of 
operation at the point of departure. 

For example, consider the following sequence: 
Location Contents 

n double precision instruction (-) 

n+1 double precison instruction (-) 
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'1+2 double precision instruction (-) 

n+3 normal 650 instruction (+) 

(the contents of n+3 are interpreted as a normal 650 instruction including the instruction 

address for sequencing. Normal execution of instructions continues until 0004 is used 

as an instruction address, at which time the program returns to the double precision 

mode beginning with the instruction in location n+4). 

Number Form 

The number form is as follows: 
. xxxxxxxxxxxxxxxxxxxx=A^ + A„. 

The low order digits of the number are stored in the location immediately following 
the one containing the high order digits. 

Instruction Form 

The instruction form is as follows: 




• double precision command (00-04) 
- 1st address A 
2nd address B 



location i 'lS'''xxxx» 



Operations 

The following is a list of the operations with their codes and estimated average 
time. 

Code Operation Estimated Average Time (ms) 

00 B=0, C — >A 23.88 

01 A+B — *.C 49. 46 

02 A-B-^C 49.46 

03 AxB-^C 113.40 

04 A/B-^C 182. 52 

Since the factors are assumed to be less than one, the numbers must be scaled 
so that the results will be less than one. If any answer should exceed one, the machine 
will stop. Also, the condition |A1 <|Blmust be satisfied for the divide operation. 

This system requires 115 storage locations. It uses locations 0000 through 0114. 
^1/ (^very location within this block is used, thus making it easy to incorporate this 

program with other programs.) The translating routine, TRl, may be used to translate 
this routine to any desired locations. It should be translated by an even amount, however, 
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to preserve the even-odd conditions. 

Explanation of Programs. 

The interpretative routine determines where the next instruction is located, obtains 
that instruction, and analyzes it to see if it must be interpreted or if it is a normal 
650 instruction. If it is a normal 650 instruction it is executed as such. If it must 
be interpreted the interpretative routine continues to analyze it, obtains the required 
factors, transfers control to the proper sub-routine. 

The constant in location 0091 facilitates use of TRl. During read in by TRl this 
constant is modified so that its instruction address digits contain the amount of 
translation. This is then added to the operation code so thai control will be 
transferred to the translated sub-routine during analysis of the operation code. 

The store, add, subtract, and multiply sub- routines are straightforward and 
simple. The latter part of the add sub-routine is common to both the subtract and 
multiply sub-routines thus saving memory locations. 

The divide sub-routines uses the following approximation to perform double 
precision division. 

A A 
"B"" B1+B2 




A l-Bp 

The sub-routine computes -^ and -^ . The multiply sub-routine is used to 

multiply these two factors to obtain the final quotient. 
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IBM IBM TYPE 650 PROGRAM SHEET 

TimlP-Miirk 

PRORI FM. INTFIRPRKTATrON WRITTEN BY:. 



FORM NO. 22-6181-0 
PniNTED IN U.S.A. 



LOCATION 

OF 

INSTRUCTION 


OPERATION 


ADDRESS 


REMARKS 


ABBRV. 


CODE 


DATA 


INSTRUCTION 


OC04 


RAL 


65 


0007 


0061 


^ 


0061 


AL 


15 


0014 


8002 


[ Increase n to n+1 and eet instruction 


8002 


SU 


11 


n+1 


0008 


J n+1 


0008 


BRMIN 


46 


0072 


0099 


Is n+1 a normal 650 instruction? 














0072 


AU 


10 


8001 


0045 


\ 


004f5 


ATI 


10 


8001 


0053 


/ n+l is a nnrmal 650 Instriir.Hon 


0053 


RT'T. 


?,n 


onn7 


«003 




8003 


Normal 


fiSflUn 


strur.tion 




/ 














0099 


STL 


20 


0007 


0010 ,, 


n+1 must be Interpreted 


0010 


AU 


10 


0013 


0017 




0017 


RAL 


65 


8003 


0025 


Is nperation code 00? 


0025 


BROV 


47 


0028 


0032 


i. 


0028 


i.n 


69 


0081 


0084 


\ 


0084 


ST DA 


22 


0087 


8001 




8001 


T.n 


fi9 


Al 


nofi? 




9961 


STn 


?4 


0023 


no2fi 




0026 


AL 


15 


0029 


0033 


Get (A) and store it in 0023. 0031 


0033 


LD 


69 


0036 


0039 




0039 


ST r>A 


?,?. 


0093 


8001 




8001 


T.n 


69 


Ao 


0071 


1 


0071 


STD 


24 


nffsi 


0034 


' 


0034 


SLT 


35 


0004 


0046 





0046 


LD 


69 


0049 


0052 




0052 


RTDA 


9.2 


nnns 


8001 




8001 


LTi 


69 


Rl 


0115 


. 


0115 


STD 


?A 


onofi 


0009 


Ofit {R) and store it in 0006, 0015 


nonfl 


AT. 


1R 


nnfi?. 


0018 




0018 


I.n 


69 


0021 


0024 




0024 


ST DA 


22 


0077 


8001 




8001 


T.n 


69 


Ro 


0110 




Olio 


STD 


24 


0015 


0068 


f 


0068 


SRT 


30 


0009. 


0075 


1 Analyze operation code and 


0075 


ATI 


10 


0091 


8003 


\ transfer to proper s?ub- routine. 


800S 


NO OP 


on 


oono 


Oppr, 


) 






































0007- 


iin; 


0008 






AddreRR nf Instruction being 1nti»rpreted 


0014r 


00 0001 


0000 








0013- 


99 nnnn 


nnnn 








0081? 


69 0000 


0067 








0029: 


00 0001 


0000 








0036- 


MJmSL 


0071 






Constants 


0049: 


69 0000 


01TS 








00R9.: 


00 0001 


onoo 








0021: 


69 0000 


0110 








0091: 


00 0000 


mmiL_ 






f 
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IBM IBM TYPE 650 PROGRAM SHEET 
PRORI FM- 00: (C) >. A WRITTEN BY: 



FORM NO. 22-6181-0 
PRrNTED IN U.S.A. 



LOCATION 
INSTRUCTIOfI 


OPERATION 


ADDRESS 


REMARKS 


ABBRV. 


CODE 


DATA 


INSTRUCTION 


0032 


AT. 


IB 


0037 


no4i 




0041 


i.n 


fifl 


0044 


8009. 




Rnn5! 


ST D 


24 


Al 


3106 


Transfer contents of C to A 


0106 


AL 


15 


0{t43 


104R 




0048 


LD 


fi9 


0051 


8002 


\ 


8002 


ST n 


714 


Ao 


)004 


) 


















































0037: 


25 0000 ( 


106 






1 Constants 


0043: 


00 0000 J 


'898 





















IBM IBM TYPE 650 PROGRAM SHEET 

Tr».lc-Mark 

PRORI FM ftl- (A)i(R) >r WRITTEN BY: 



FORM NO. 22-8181-0 
PRINTED IN U.S.A. 



LOCATION 



sTijr 



INSTRUCTION 



OPERATION 



ABBRV. 



CODE 



ADDRESS 



DATA 



INSTRUCTION 



REMARKS 



-mm- 



Tonr 



0006 



0011 



0011 



AL 



15_ 



nnifi 



noi9 



0019 



■Mi- 



iO- 



0023 



sioa^ 



. Aj + Bi + A2 + R2 



QQ27 



AL_ 



IS 



0031 



0035 



0035 



BRQV 



41. 



0107 



0040 



TftRt fnr nvftrflnw 



0040 



STU 



21. 



0044 



0047 



Store siiTT) In C 



04 7 



5ILi_ 



iQ_ 



0051 



imQ4_ 



0107 



STOP 



Sil. 



0000 



0107 



Overflow nn addition (or subtraction) 



IBM IBM TYPE 650 PROGRAM SHEET 

Tii.ilc-.Mn(k 

PRORI FM. 02; (A)-(B)^C WRITTEN BY:. 



FORM NO. 22-6IBI-0 
PRINTED IN USA. 



LOCATION 
INSTRUCTION 


OPERATION 


ADDRESS 


REMARKS 


ABBRV. 


CODE 


DATA 


NSTRUCTION 


0001 


RSII 


fil 


0006 


0012 


\ 


0012 


SL 


16 


0015 


0010 
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A METHOD FOR PERFORMING COMPLEX ARITHMETIC ON THE IBM TYPE 650 

G. R. Trimble, Jr. 

Introduction 



This method for performing complex arithmetic on the 650 has been developed 
since it is sometimes necessary to perform computations using complex numbers. 
The basic operations of addition, subtraction, multiplication and division using 
complex numbers are performed by this method. It was designed with coding 
convenience as a prime objective and uses an interpretative routine to perform 
complex arithmetic. 

Complex instructions, that is, instructions which are to be interpreted, have a 
negative sign. A complex instruction consists of a 2 digit operation code, a 4 digit 
address specifying the location of the first factor and another 4 digit address 
specifying the location of the second factor. The result of the operation is always 
stored in a pair of memory locations referred to as C. A 5th operation, "Store", 
is used to transfer the result from C to the locations desired by the programmer. 
In order to store the result where desired, it will usually be necessary to follow one 
of the 4 basic arithmetic operations by a store instruction. The net effect of such 
a procedure is four 3-address operations. 

The address given refers to the real part of the factor to be operated upon. The 
interpretative routine automatically selects the imaginary part from the next memory 
location. Thus, the real and imaginary parts of a complex number will always be 
stored in successive memory locations. 

Should a positive instruction appear, it is interpreted as a normal 650 instruction 
and subsequent instructions are not interpreted. Thus, when a positive instruction 
occurs the 650 operates in its usual D-I mode. This continues until an instruction 
address of 0054 is given which causes the control to return to the complex mode of 
operation at the point of departure. 

For example, consider the following sequence: 
Location Contents 



ii complex instruction (-) 

"^+1 complex instruction (-) 
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n+2 complex instruction (-) 

n+3 normal 650 instruction (+) 

( the contents of n+3 are interpreted as a normal 650 instruction including the instruction 

address for sequencing. Normal execution of instructions continues until 0054 is used 

as an instruction address, at which time the program returns to the complex mode 

beginning with the instruction in location n+4). 

Number Form 

The number form is as follows: 

. XXXXXXXXXX+. xxxxxxxxxxi=A J + A-l 

The imaginary part of a number is stored in the location immediately following 
the one containing the real part. 

Instruction Porm 



The instruction form is as follows: 



' complex command (00-04) 
■ 1st address A 
-2nd address B 



location i 'STtxxxk'icxxx^ 

Operations 

The following is a list of the operations with their codes and estimated average 
time. 

Code Operation Estimated Average Time (ms) 

00 B=0, G^vA 23.50 

01 A+B >C 45.64 

02 A-B^ yc 45.64 

03 AxB >.C 96. 50 

04 A/B ^C 175.42 

Since the factors are assumed to be less than one, the numbers must be scaled 
so that the results will be less than one. If any answer should exceed one, the machine 
will stop. Also, the condition |A| <|B1 must be satisfied for the divide operation. 

This system requires 126 storage locations. It uses locations 0000 through 0127. 
Every location within this block is used, thus making it easy to incorporate this 
program with other programs. The translation routine, TRl, may be used to translate 
this routine to any desired memory location. It should be translated by an even amount 
however, to preserve the even-odd conditions. 
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Explanation of Program 

The interpretative routine determines where the next instruction is located, obtains 
the instruction and analyzes it to see if it must be interpreted or if it is a normal 650 
instruction. If it is a normal 650 instruction it is executed as such. U it must be 
interpreted the interpretative routine continues analyzing it, obtains the required 
factors, and transfers control to the proper sub-routine. 

The constant in location 0091 facilitates use of TRl. During read in by TRl this 
constant is modified so that its instruction address digits contain the amount of 
translation. This is then added to the operation code so that the control will be 
transferred to the translated sub-routine during analysis of the operation code. 

The sub-routines are straightforward and simple. The latter part of the add 
sub-routine is also used by the subtract sub-routine to save memory locations. The 
divide sub-routine may require that the dividend and divisor both be shifted right one 
place if Bi + Bl should exceed one. 
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